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Nuclear gluon densities are of great importance to the physics of relativistic heavy 
ion collisions, in particular, in assessing the origin of J/^-suppression. We describe our 
attempts to distinguish various models of the gluonic EMC-effect, using the existing 
J/^-production data in proton-nucleus collisions. We find that no model is capable of 
explaining all the features of the high precision E772 data although the overall trend 
suggests this to be more a matter of fine-tuning the models than the presence of new 
physical effects. 
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1. INTRODUCTION 



Theoretical investigations of several aspects of the heavy-ion collisions at high energies 
depend crucially on the precise knowledge of the gluon distributions of the colliding nuclei. 
A well-known prime example is the //-^-suppression signal []T],|2| of quark-gluon plasma 
(QGP). It has been shown|3]||] that this signal and its transverse momentum dependence 
can be mimicked by the nuclear structure function effects. The nuclear glue can thus be 
a potentially serious source of background which has to be subtracted off to see QGP. 
Another important reason to size up the nuclear glue is the EMC-effect observed in the 
ratio of quark structure functions, p = (x) / 'AF2(x) . There are a variety of diverse 
theoretical ideas which "explain" the data on p. Since most of these result in a prediction 
for the corresponding gluon ratio, p g , any attempt to obtain experimental information on 
p g is likely to assist in distinguishing between the various models and in pinning down 
the likely cause of this effect. 

Extraction of gluon densities from the data is known to be a difficult exercise al- 
though there are various known methods such as the direct photon production or the 
//^-production. The problem is that in most cases the quark distributions also con- 
tribute and one has to invent clever tricks or resort to some approximations to separate 
the gluon densities. Recently, Drees and Kim advocated the study of associated pro- 
duction of J / ip and a photon at large px as a clean probe of the gluon density in pp and 
ep collisions. We extended || their considerations to pA and A A collisions for the fixed 
target experiments at FNAL and the collider experiments at RHIC. We found that the 
rates for these processes at both FNAL and RHIC are substantial and a large range of 
x is covered by combining them, thus enabling a cleaner and better determination of p g . 
For more details, including a technique to improve the statistics, we refer the reader to 
Ref.§. 

2. J/^-PRODUCTION 

While the associated production of J/ip and photon at large px could serve as a good 
probe of p g in future, one naturally asks the question whether the existing data can in 
some way be used for this purpose. In particular, the recent high statistics J/ip and 
T- production data from the E772 experiment^, combined with the older data || have 
been exploited by Gupta and Satz|| to obtain p g . We argue that their assumptions are 
not valid and one therefore needs to take a less ambitious approach of 

checking whether the data can be explained by the existing models for p g . Before we 
do so, let us recall briefly the theoretical formalism. 

2.1. Theoretical Models 

The basic perturbative QCD-subprocesses which contribute to the production of J/ip 
are 

99 -> cc; qq -> cc and gg -> ccg; qq -> ccg; qg -> ccq; qg -> ccq (1) 



The first two of these contribute at 0{a 2 s ) while the rest do so at O(a^). Furthermore, 
since the cc pair will hadronize to J/ip, it is also clear that the former will result in J/ip 
at low (~ 0) px while the latter will yield large px for J/ip. Two popular [|l(J models for 



hadronization of the cc pair to J/ip are 1) semilocal duality and 2) colour singlet model. 
The former simply consists of computing the cc pair cross section with the appropriate 
perturbative QCD cross sections and the chosen quark and gluon densities and multiplying 
them by a constant if the pair mass lies in the range between Mj/^ and Mp>- Hadronization 
in the latter model, on the other hand, is governed by the quarkonium wavefunction (or 
its derivative) at the origin. An intersting feature to note is that the charmonium is 
produced at the basic QCD vertex in this model and the kinematics is very different. 
Thus, e.g., the J/ip at large pp is produced in a basic 2^2 process in contrast to the 
semilocal duality model above. 

2.2. Results 

In order to extract p g , Ref. assumed that the data can be described by 0(a 2 s ) 
partonic cross section and by further assuming that gluons dominate. Since the pr-range 
of E772 goes up to 2.25 GeV, the first assumption implies rather high values for primordial 
Pt'- (pt) ~ 1 — 2 GeV. Furthermore, as Fig. [1] shows, the assumption of gluon dominance is 
both x^-dependent and structure function parametrization-dependent. What is plotted in 
the figure is the ratio of the product of quark structure functions to that of gluon structure 
functions which act respectively as the weights of the quark-antiquark and gluon-gluon 
terms in pp collisions at r = Mj/^/s/s = 0.0775. The various curves are for different 
popular parametrizations of the structure functions. The corresponding curves for T cast 
even more doubt on the validity of gluon dominance. 

We therefore used the 0(a^) partonic subprocesses to compute the pr-distribution of 
J/ip and T for the E772 acceptance for both the semilocal duality model and the colour 
singlet model. The details of our calculations can be found in Ref. |TTJ. Here we wish to 



make only some general comments. Since the E772 data is unfortunately not sufficiently 
differential, we had to extend our integration range over larger xp where our approach 
may be inadequate. Fortunately, the contribution from the large xp region is small and 
thus will not affect our conclusions. Since the rcp-spectrum of E772 also includes low p T 
data, we are unable to compute it at all in our approach. Figs. |^ and ^| exhibit our results 
along with the E772 data for their lightest and heaviest target for three different models 
of p g . Again, more details about these models, including original references can be found 
in Ref. [IT]. What is remarkable is that for Tungsten, none of the model agrees with the 



data, although for carbon, the so-called Gas model seems to perform well. For the other 
two targets of E772, situation is akin to that in Fig. || 

Performing the same computation for the colour singlet model is computationally a lot 
easier since the final state is simpler. In order to highlight the independence of the above 
results on the hadronization models in a typical example, we display in Fig. |] the ratio of 
the results for the two models discussed above as a function of pp for the Tungsten nucleus. 
The agreement is really remarkable, especially when one takes into account the differences 
in kinematics as well. We conclude from here that the data are really a reflection of the 
structure function effects alone, although the models of p g considered here need to be 
fine-tuned. 

Essentially the same conclusion seems to emerge from our analysis of the E772 data on 
T-production. Since only a(pp) as a function of pp is available in this case, where a{pp) 
is obtained by fitting the cross section at each value of pp for all nuclei to the form A a , 



we followed the same E772 procedure and obtained a(pr) in the range 1-4 GeV. As seen 
in Fig. a very good agreement with the E772 data is obtained for both the Gas model 
and a 6-quark cluster model for pt upto 3 GeV. Surprisingly, the computations remain 
much flatter at even higher pt values whereas the data show a dramatic increase. Since 
the large p r -domain is theoretically better suited for the description employed above, this 
discrepancy is even more striking and needs to be understood more carefully. 
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Figure Captions 

Figure 1. Test of gluon dominance as a functionofxF. 



Figure 2. R(pr) vs. pr: Carbon 



Figure 3. R(pt) vs. p^: Tungsten 



Figure 4. Ratio of models: Tungsten 



Figure 5. a(pr) vs. px- T 
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